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ABSTRACT 



Context. The Kepler space mission is reaching continuous observing times long enough to also start studying the fine structure of the 
observed pressure-mode spectra. 
<~i ' Aims. In this paper, we aim to study the signature of stellar evolution on the radial and pressure-dominated / = 2 modes in an ensemble 

CL(, of red giants that show solar-type oscillations. 

Methods. We use established methods to automatically identify the mode degree of / = and 2 modes and measure the large (Avt.) 
and small {Svqi) frequency separation around the central radial mode. We then determine the phase shift e^ of the central radial mode, 
i.e. the linear offset in the asymptotic fit to the acoustic modes. Furthermore we measure the individual frequencies of radial modes 
^ \ and investigate their average curvature. 

Results. We find that e^ is significantly different for red giants at a given Av^ but which bum only H in a shell (RGB) than those that 
have already ignited core He burning. Even though not directly probing the stellar core the pair of local seismic observables (Avc, e^) 
can be used as an evolutionary stage discriminator that turned out to be as reliable as the period spacing of the mixed dipole modes. We 
^ [ find a tight correlation between e^ and Av^ for RGB stars and unlike less evolved stars we find no indication that e^ depends on other 

^^ , properties of the star It appears that the difference in e^. between the two populations becomes smaller and eventually indistinguishable 

ffN ■ if we use an average of several radial orders, instead of a local, i.e. only around the central radial mode, large separation to determine 

the phase shift. This indicates that the information on the evolutionary stage is encoded locally, more precisely in the shape of the radial 
rf^ [ mode sequence. This shape turns out to be approximately symmetric around the central radial mode for RGB stars but asymmetric 

• , for core He burning stars. We computed radial mode frequencies for a sequence of red-giant models and find them to qualitatively 

C^ ■ confiiTn our findings. We also find that, at least in our models, the local Av is an at least as good and mostly better proxy for both the 

'<— ^ ' asymptotic spacing and the large separation scaled from the model density than the average Ay. Finally, we investigate the signature 

Cn ' of the evolutionary stage on Svoi and quantify the mass dependency of this seismic parameter 

T-H ' 

Key words, stars: late-type - stars: oscillations - stars: fundamental parameters - instrumentation: Kepler 
^ , 

X 

5— I 1 . Introduction to expand. Because of the drop in surface opacity and the high 

. ?.' sub-surface opacity, the subsurface layers become convective in 

Once a main-sequence star has depleted its supply of H in the o^der to maintain the flow of energy out of the star. The star's 15 

core it begins to bum H in a sheU surrounding the nearly isother- luminosity and radius increase while its suiface temperature re- 

mal He core. With the mean molecular weight increasing, the ^ains relatively fixed, i.e., the star climbs the RGB. Eariy in this 

5 core is forced to contract under gravity and heats up. To maintain ph^se the convective envelope extends deep into the interior but 

the increasing flow of power from the central regions, the outer ^s the star evolves, the convective region retreats outwards, 
layers expand and cool down until the surface opacity dramati- Eventually the core is hot enough to ignite He (triple-alpha) 20 

cally drops due to the shortage of available free electrons from ^^^. ^^^ ^^^ ^^^^^ contracts. In lower mass stars, the cen- 

lonised inetals that contribute to the H" opacity. The star has ^^^j ^^^^-^-^^ ^^ ^-^ ^^^^^ ^^^^ ^^^ ^^^^^^^ -^ degenerate. 

1 reached the red giant branch (RGB) phase. The core continues to ^j^^ ^^^^ temperature dependence of the supporting electron gas 
contract, its temperature continues to nse, the luminosity of the ^.^ j^^^^ ^^ ^ ^^^^.^^^ runaway, with the He burning rate in- 

burnmg H shell continues to rise, and the outer layers continue ^^^^^-^^ ^^^-^^^^ unconstrained by pressure expansion until de- 25 

generacy is removed. The star quickly settles onto the horizontal 

Send offprint requests to: thomas.kallinger@sterkuleuven.be branch (also called the red clump, RC), burning He in the core 
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and H in a shell, with the latter providing most of the energy. 
Eventually, the star will exhaust He in the core and the structural 
30 changes that occurred earlier will repeat, this time entering the 
asymptotic giant branch (AGB) phase, with the star moving back 
up the giant branch. The star is then burning He in a shell sur- 
rounding an inert carbon core, and burning H in a shell slightly 
farther out, being still the dominant energy source. 

35 Even though the general concep t of red-giant evolu tion is 

well understood (e.g., Iben 1968; Swe igart & GrossI 1978b . there 
are a number of outstanding problems, such as mass-loss on the 
giant branch, deficits in modelling the convective energy trans- 
port, and the internal mixing from rotation, that limit our un- 

40 derstanding of the chemical evol ution of galaxies, w ith conse- 
quences even for cosmology (e.g. JChristensen-Dalsg aard 2009]). 

A problem in addressing these questions is our inability 
to distinguish the different evolutionary stages of red giants. 
Having quite different internal structures, they still cover the 

45 same region in the Hertzsprung-Russell (HR) diagram and we 
cannot distinguish them based on their surface properties, i.e., 
with classical methods like spectroscopy. Asteroseismology of- 
fers a solution, being an elegant tool to probe the internal struc- 
ture of a star by making use of its natural eigenfrequencies, 

50 which depend on the physical properties of its interior (e.g., 
lAerts et al...2010i) . 

All cool stars with a convective surface layer are believed 
to show high-overtone oscillations of low spherical degree I. 
These so-called solar-type oscillations are intrinsically damped 

55 and stochastically excited by the turbulent flux of the near- 
surface convection. For main-sequence stars, they are pressure 
(p) modes that are mainly confined to the outer envelope and 
which follow a comb-like pattern in the frequency spectrum. The 
existence of solar-type oscillations in the more evolved red gi- 

60 ants is n ow well established observationa lly, most recently f rom 
CoRoT (jBaglin et al.ll2006h and Kepler OBoruckiet al.ll2010l) . 

However, the dense He cores of red giants complicate the 
situation, allowing the existence of gravity (g) modes at fre- 
quencies that overlap the p-mode frequency spectrum. Red-giant 

65 models (e.g., Guenther et al. 2000; Dziembowski et al. 2001) 
predict a very dense spectrum of non-radial mixed modes with 
a g-mode character in the core and a p-mode character in the 
envelope, whose regular period spacing is directly dependent on 
the b uoyancy frequency in the core region (e.g., IDziembowskil 

70 11977) . However, g-dominated mixed modes have very high in- 
ertias, leading to very small, often unobservable, amplitudes at 
the surface. From a theoretical perspective, however, some of 
the mixed modes have substantial amplitudes in the envelope 
due t o resonant co upling between t he g- and p- mode cavities 

75 (e.g., IChristensen-Dalsgaard 2004: Dupret et all 12009) . Their 
mode energy is then concentrated in the envelope giving rise 
to th e surface amp litudes and making them clearly detectable 
(e.g.,lKallinger et al. 2008; De Ridder et al. 2009; Bedding et al. 
12010 ). For modes with / > 1, mostly single modes are observ- 

80 able per radial p-mode order with their frequencies following the 
radial mode pattern. Only for dipole modes (1=1) the trapping 
is less selective and more than one mode reaches sufficient sur- 
face amplitudes. They appear in groups in bet ween consecutive 
radial modes forming a characteristic pattern (iBecket al.ll201lh 

85 with an average period spacing that is smaller than those of the 
g-dominated mixed modes that are confined to the core. Thanks 
to the quality of the data obtained with NASA's Kepler space 
telescope for h undreds of red giants . Bedding et al. (201 1 ) (and 
simultaneously iMosser et a0201 lab found that the observed av- 

90 erage period spacings are still a valuable probe for the core prop- 



erties. They permit a distinction between evolutionary stages of 
otherwise indistinguishable red giants. 

Here, we concentrate on radial modes and how they are af- 
fected by the evolutionary stage of the star Their frequency pat- 
tern is largely defined by a characteristi c spacing and a phase 95 
shift. From 5 -month long CoRoT data, iMosser et aP (1201 Ibl) 
postulated an universal pattern for red-giant oscillations. They 
measured an average large frequency separation (Av) and in- 
dicated that the absolute p-mode frequencies of all red giants 
are basically a function of Av and that there is no difference 100 
in the p-mode spectrum of a RC star and an RGB star if they 
have the same Av. On the other hand, iBeddin^ et al.l (1201 ih 
found evidence that the frequency spectrum of a red giant with 
a given Av is shifted in frequency according to its core proper- 
ties, i.e. the phase shift of an RGB star differs from that of a RC 1 05 
star, clearly contradicting the universal pattern. We find that this 
contradiction is caus ed by different definitions of Av. Whereas 
IMosser et al.l (l2011bh determined an average Av of all observed 
p modes and took in to account a gradient of Av with frequency, 
iBedding et al.l (1201 1*) adopted the large frequency separation that 1 1 
best aligns the strongest radial modes in an echelle diagram, 
which is often different from an average value of all modes. 
Both methods are potentially sensitive to observational charac- 
teristics, such as the signal-to-noise ratio, making it difficult to 
use them for a high-precision comparison of different stars. We 1 1 5 
re-examined the Kepler data using a more distinctive "local" def- 
inition of Av, namely the average separation of the three central 
radial modes. Indeed, if we compare stars with the same cen- 
tral Av, we find that the central radial modes of core He-burning 
stars have a smaller phase shift than do stars still burning only H 1 20 
in a shell. Therefore, the three central radial modes can be used 
to distinguish between the different evolutionary stages of red 
giants. 

We also investigated the fine structure of the p-mode spec- 
trum and found that the radial mode sequence of RGB stars is 1 25 
approximately symmetric around the central radial mode in con- 
trast to those of RC stars, which is asymmetric or even antisym- 
metric. Also the small frequency separation between adjacent / 
= and 2 modes is affected by a star's evolutionary status, lead- 
ing to a small frequency separation th at is on average large r for 1 30 
RC stars than for RGB stars (see also ICorsaro et al.ll20T2l) . We 
checked the seismic parameters for a scaling with stellar mass 
and found it to clearly affect the small spacing. Finally, we com- 
puted radial mode frequencies for a sequence of red-giant mod- 
els and found them to qualitatively confirm our findings. 135 

2. Observations 



The Kepler space telescope (iBorucki et al.l 1201 Ot iKoch et aP 
2010) was launched in March 2009 to search for transiting Earth- 
sized planets in and near the habitable zones of Sun-like stars. 
Kepler houses a 95-cm aperture modified Schmidt telescope that 1 40 
points at a single field in the constellations of Cygnus and Lyra, 
feeding a photometer with a wide field of view to continuously 
monitor the brightnesses of about 150000 stars. This makes it 
also an ideal instrument for asteroseismology and the Kepler 
Asteros eismic Science Conso rtiunJ] has been set up to carry this 1 45 
out ("see lGilUland et al]l2010l for an overview and first results). 

Kepler observations are subdivided into quarters, starting 
with the initial commissioning run (lOd, QO), followed by a 
short first quarter (34 d, Ql) and subsequent full quarters of 
90 d length. Our studies are based on the long-cadence (29.42 150 



http://kepler.asteroseismology.com 
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min sampling: lJenkins et al.ll2010h data spanning from QO to Q7 
with a total of about 27 200 measurements. Apart from occa- 
sional losses of the satellite's fine pointing and scheduled re- 
orientations of the spacecraft, the ~599 days-long observations 
were continuous, with in an overall duty cycle of about 93 %. 
The Kepler raw da ta were reduced in the manner described by 
iGarcia et alJ (12011 ^) and subsequently smoothed with a triangu- 
lar filter to suppress residual instrumental long-term trends with 
time scales longer than about 10 d. The power density spec- 
tra were normalised to recover the full sine-amplitude of an 
injected signal. Our sample includes a total of 923 red giants 
showing a clear power excess due to pulsation and has already 
been used for other studies, such as the inves tigation of rela- 
tions between the d iff'erent seismic (Huberet al. 2010t 201l1: 
lMosseretaill2012al) and granulation (iMathuret al.ll201 1l) ob- 
servables, the seismic d etermination of fundamental parameters 
dKallinger et a l. 2010b, hereafter referred to as Paper I), a com- 
parison o f global oscillation p arameters derived from different 
methods (iHekker et al.l 201 \m. and the analysis of non-r adial 
mixed modes dBeddinget al J 120111: iMosser et all (20 12b') and 
how t hey are used to constrain the core-rotation rate (3eck et al., 
I2012h . 



3. The acoustic p-mode spectrum 

The observed oscillation modes of solar-type pulsators are typi- 
cally high-order acoustic modes. If interaction with g modes can 
be neglected, the linear, adiabatic, high-order acoustic modes in 
a spherically symmetric star follow an asymptotic relation for 
the frequencies with a radial ord er n and a spherical degree / 
(I Vandakurovl 1 9671: iTassouJI 1 9801) : 



v„,/ ^ A I n + - + e I - (io/. 



(1) 



where the characteristic (or asymptotic) spacing^ A ^ v„+i,; - v„j 
is the inverse sound travel time across the stellar diameter, es- 
senti ally scaling with the squar e root of the mean stellar density 
(e.g.. iBrown & Gillilandlfl994l) . The phase shift e is a dimen- 
sionless ofi^set and doi denotes the small frequency separations 
of non-radial modes relative to radial modes, which is given by 
v«,o - >'n-i,2 in the case of / = 2 modes. In this first-order ap- 
proximation, the modes are expected to fall al ong vertical ridge s 
when folded with A in an echelle diagram (iGrec et al.llI983l) . 
with the radial modes shifted by e times A from the origin and 
the non-radial modes shifted by d()i relative to the radial modes. 

185 This simple concept appears to be quite promising because it 
suggests that some observed frequency separation can directly 
be related to physical quantities of a star. However, in reality the 
situation is more complex. The actual modes fall along curved 
ridges in the echelle diagram and the departure from the asymp- 

190 totic relation is afunctionoffrequency(orequivalentlyofn) pre- 
dominantly determined by the mean stellar structure and atmo- 
sphere. Additionally, sharp changes in the acoustic variables, e.g. 
at the boundary of the He II ionisation zone, cause small quasi- 
periodic variations aroun d t he curved mode seq uences (see, e.g. 

195 iHoudek & Goughll2007l and lMiglio etani201ol for the Sun and 
a red giant, respectively, and iHekker et al.l 1201 lal for red-giant 
models). 

In the classical asymptotic theory, the higher-order correc- 
tion terms are presented as functions of v. The radial mode pat- 



tern can be written as v(n) - A(v)[n -i- e] with e fixed, or math- 
ematically equivalent as v(n) = A[n + e(v)] with A fixed and 
preferably e(v) - ef{v). We prefer the latter because it allows A 
to have a single value th at is independent of both / and frequency 
(see also iBedding||20I 11) . This can be developed as follows. The 
frequencies have to satisfy t he Eigenfrequency equation (e.g., 
[Roxburgh & Vorontsovll2000l) : 



^«,/T = ^ n + 2 I + "''(^^ ~ ^'(■^)' 



(2) 



where ai(v) is an outer phase shift primarily determined by the 
mode propagation in the outer layers (including the atmosphere) 
of the star, and 6i(v) an inner phase shift determined primarily 
by the structure of the inner regions. The Eigenfrequency equa- 
tion can be rewritten so that the phase shift and small separa- 
tion from Eq.[T]are expressed in terms of phase shift differences 
(Roxburgh, private communication): 



e(v) = -\aoiv) - doiv) 



doi(v) = -| Q'o(v) - ai(v) + 6i{v) - 6q(v) 



(3) 



(4) 



^ Note that the asymptotic spacing is usually referred to as Av but in 
order to keep it clearly distinct from any observed quantity we prefer A 



In fact, the outer phase shifts ai(v) are almost equal for differ- 
ent / and can be taken as independent of I. Consequently, d()i{v) 
depends only on the structure of the inner region, more explic- 200 
itly on the sound-speed gradient near the stellar core. This arises 
because the small spacing takes the frequency difference of two 
modes with nearly identical eigenfunction shapes in the outer 
layers of the star. Therefore, it is sensitive to the structure of 
the deep interior The phase shift e(y) has contributions from the 205 
outer and inner regions. For the Sun, a'o(v) is ~5 times larger than 
6o(v) (Roxburgh 2010) so that the dominant contribution for the 
solar modes comes from the surface. This is why e is sometimes 
referred to as the "surface term", which we find not quite appro- 
priate because it clearly has also contributions from the interior 210 
Obviously, the departure from the first order asymptotic rel ation 
is given by the variation of oq - 6o with frequency. 

A problem is to relate the observed large (Av) and small 
(6vq2) frequency separations and phase shift (e) to the coeffi- 
cients discussed above, which are not straight-forwardly accessi- 21 5 
ble. The large frequency separation, for example, should ideally 
be measured from the modes closest to the acoustic cutoff fre- 
quency, which presumably follow the characteristic spacing best. 
But these modes often have very small amplitudes or are unob- 
servable. An average Ava of all observable (radial) modes might 220 
be a good alternative (e.g.. White et al. 2011b). Even though it 
can be determined more accurately and should average out the 
high-frequency variation along the mode sequence it is poten- 
tially sensitive to the mode curvature. For the Sun, for example, 
the separation of consecutive radial modes varie s by almost 6 % 225 
over the observed range (e.g.. lChapIin et al. 1999). with the local 
separation dropping towards higher frequencies for the lowest 
radial orders, staying relatively constant around the centre of the 
power excess hump and slowly increasing towards the highest 
radial orders. Hence depending on where exactly Ava is mea- 230 
sured in this sequence one might end up with values that are sig- 
nificantly different. A careful inspection of many red-giant fre- 
quency spectra has shown that the curvature of th e radial mode 
sequence is small (on the order of 1 % of Av, e.g..' Mosseret al.l 
2011b) but can be quite different in amount or even direction. 235 
More interestingly, it appears that the curvature always has a 
turning point that is near the centre of the power excess, i.e. the 
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curvature is "in phase" at this distinctive point. This is also the 
case for the Sun and other solar-type oscillators. 

To see how well the different choices of Av represent the 
asymptotic spacing we investigate theoretical frequencies for a 
sample of representative models in Sect. l4.1l 

3. 1 . Practical implications for red giants 



295 



iMosser et al.l ( 1201 Ibl) found empirical relations between e, 6v, 
and Av for the red giants observed by CoRoT. Based on these 
scalings they constructed theoretical oscillation patterns and cor- 
related them to the observed one to precisely refine their estimate 
of an average Av comprising all (except dipole) observed modes. 
In their approach they assumed that the absolute position of the 
radial-mode pattern depends only on Av, i.e. with e being a func- 
tion of only Av, regardless of whether the star burns He in the 
core or still only H in a shell. Consequently, all red giant oscil- 
lations should follow a universal pattern and we should not see 
any difference in the radial mode spectrum of RGB and RC stars 
if they have the same Av. This can be seen in Fig.[Tl where we 
plot normahsed echelle diagrams of 923 red giants observed by 
Kepler folded with an average Av measured using the method of 
iMosser et al.l (l2011 b) and stacked on top of each other with the 
resulting total power colour-coded. All spectra show the same 
pattern with the / = and 2 (and 3) modes falling on top of each 
other and the phase shift following a simple trend. We note that 
the scatter along these ridges does not only reflect observational 
uncertainties or the limited lifetimes of the modes but also the 
fact that more than one radial order per stars contributes to the 
total power As expected the / = 1 ridge is somewhat broadened 
due to multiple peaks per radial order 

The situation changes if we use a more local Av instead of an 
average one to fold the spectra. This can be seen in the left panel 
of Fig. 121 Whereas the observed modes of stars with their oscilla- 
tion frequencies above about 50yuHz form a regular pattern with 
the radia l modes (righ tmost ridge) following the empirical re- 
lation of IMosser et al.l (i2011b) (white dashed line), stars below 
about SOyuHz show a less clear structure (see also iHuber et all 
12010 *). The latter is the range where we expect both RGB and RC 
stars (e.g. Miglio et al. 2009). To check if this structural break is 
due to different evolutionary stages of stars contrib uting to this 
diagram, we adopted the identification of Bedding et al.l (1201 ih 
and indeed found that if we plot only RGB stars (middle panel), 
the modes form a similar pattern to Fig.[T] even in the range that 
is also covered by RC stars (right panel). We note that for RGB 
stars, modes below about 50 /kHz are barely visible in the left 
panel because this range is dominated by RC stars, which clearly 
outnumber the RGB stars. 

The different positions of the mode ridges in Fig. |2] indicate 
that RC stars have a phase shift that differs syste matically from 
that o f RGB stars. This was already noted by iBedding et al.l 
(1201 lb but their measurements of e were not accurate enough 
to serve as an unambiguous indicator for the stars' evolutionary 
stage. This might be caused by their specific choice to determine 
Av over all observed modes. 



3.2. The central large frequency separation 

To follow up on this interesting difference we consider a more 
distinctive definition of the large frequency separation, namely 
the average separation of the three radial modes around the 
centre of the pulsation power excess (Avc), as introduced by 
iKalhnger et"^ (l2010ch and Paper I. 
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Fig. 1. Stacked and normalised echelle diagrams of 923 red gi- 
ants observe d by Kepler folded w ith an average Av computed 
according to IMosser et al.l (1201 Ibl) with the total power colour- 
coded. The individual spectra are corrected for the granulation 
signal. Numbers indicate the mode degree. 



First we fitted a global model to the observed power den- 
sity spectra to determine the granulation background signal. 
The model consists of a superposition of white noise, the sum 
of super-LorentziaqJ functions, and a Gaussian to approximate 300 
the power excess due to pulsation. The latter component pro- 
vides a good estimate for the centre of the pulsation power ex- 
cess, the so-called frequency of maximum oscillation power or 
^max- For more details about the fitting we refer to Paper I and 
Mathuretal. (201ll). 305 

In the next step, we used the global model fit without the 
Gaussian component and fitted on top of that a sequence of 
Lorentzian profiles. The sequence covers three radial orders of / 
- and 2 modes around v^ax, where the mode frequencies are 
parameterised by the frequency of the central radial mode (vc o) 310 
and a large (Av^) and small (dvoi) spacing. A single linewidth 
parameter was assumed for the six profiles, which might not be 
true in reality, but has no impact on the determination of the 
spacings and significantly stabilises the fit. To account for the 
power of the intervening p-dominated mixed dipole modes, two 315 
additional profiles were used, where the linewidth was fitted in- 
dependently of the linewidth of the / = and 2 modes. Fig.|3] 
illustrates the power density spectra of RGB (top) and RC (bot- 
tom) stars with almost identical Avc, along with the best fitting 
models (red lines). In the case of the RGB star with very closely 320 
spaced dipole modes (that are almost unresolved), the two addi- 
tional profiles give a good estimate of where to expect the pure I 
= 1 p modes, parameterised by 5v()\ , which i s the separation from 
the midpoint of adjacent radial modes (e.g. Bedding et al.l2010l) . 
For the RC star, whose individual dipole modes are well sepa- 325 
rated, the two additional profiles give no useful "physical" infor- 
mation. Their purpose here is to represent the integrated power 
of the dipole modes to stabilise the fit for the I - and 2 modes. 



^ The adopted background model is a modified Lorentzian with the 
exponent equal to 4 (instead of 2) and is therefore called "super- 
Lorentzian" 



Kallinger et al.: Evolutionary influences on the structure of red-giant p-mode spectra 



200 



N 



100 



50 



all 


1 1 1 I 1 1 1 < 1 1 1 . 1 ^ 




' if . 


i 


/ 




0.2 0,4 0,6 0,. 



1 0.2 0.4 0,6 0,8 
v/hv^_ modulo 1 



1 



0,2 0,4 0,6 



0,8 



Fig. 2. Same as Fig.[T]but now folded with the central Av. The left, middle, and right panel show all, only RGB, and only RC stars, 
respectively, with the identi fication based on the g-mode period spacings. The white dashed lines represents the empirical relation 
between e and Av from Mosser et alJ (1201 Ibl) . 
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Fig. 3. Power density spectra of an RGB (top) and RC (bottom) 
star with almost identical Avc but with the central radial modes 
(Vc,o, indicated by blue arrows) at significantly different frequen- 
cies (resulting in different ff). The red lines indicate the best fit 
model to determine the p-mode spacings, Avc and 6vw.. Numbers 
give the mode degree. 



To further reduce the number of free parameters for the fit 
we parameterised the heights of the individual profiles with a 
Gaussian envelope for all modes of a given degree. The result- 
ing three Gaussians utilise a single central frequency and width 
parameter but individual heights for each mode degree. This re- 
sulted in a total of 11 free parameters to fit (compared to 24 
if everything is left free). Instead of the Markov-Chain Monte 
Carlo (MCMC) algorithm used in Paper I, we applied a Bayesian 
inference tool called MulHNest (Feroz et al. 2009). The big ad- 
vantage of the nested sampling algorithm is that it provides the 
posterior distributions for the parameter estimation as well as 



the model evidence, without the need for multiple runs (as it 340 
would be necessary for a MCMC algorithm) and thereby sig- 
nificantly outperforms any other fitting method in terms of ef- 
ficiency and robustness. With a reliable statistical measurement 
in hand, telling us how well a given model represents the obser- 
vations, we could investigate different parameterisations of the 345 
eight Lorentzian profiles. Usually, more complex models (i.e., 
with more free parameters) tend to fit some data better than less 
complex models but in a Bayesian analysis a model is assigned 
a penalty for its complexity and needs to fit the data consider- 
ably better to get a higher model evidence than a less complex 350 
model. Fitting different complex models (e.g., with independent 
mode height and/or linewidth parameters) to a number of typi- 
cal red-giant spectra, we found that the above model represents 
the observations best. Adding further parameters (i.e., making 
the model more complex) did not improve the fit (i.e., the model 355 
evidence) for the present observations. 

During the fit, the centre and width of the mode height en- 
velopes were allowed to vary over the ranges v,nax ± o"? and 
CTg + O.ScTg, respectively, where o-g corresponds to the width 
of the power excess, as measured from the background fit. The 360 
linewidth parameters were sampled in a range that corresponds 
to a mode lifetime of 10 to 100 days for the / = and 2 modes, 
and 2 to 100 days for the / = 1 modes. Note that a lifetime of 
100 days corresponds to a linewidth that is about twice the fre- 
quency resolution of the observations. The individual envelope 365 
heights were sampled between and twice the power density of 
the highest peak in the frequency range of interest. A critical pa- 
rameter is the central frequency because the central mode must 
be a radial mode in order for the fitted spacings to be meaningful. 
We allowed Vcfi to vary in the range Vmax ± 0.55 times an initial 370 
guess for Av taken from previous analysis, as in Paper I. Since 
we already have good estimates for Av and 5v()2 from previous 
work, we can now be more restrictive for their sampled parame- 
ter ranges than in Paper I. The spacing parameters were allowed 
to vary between 0.8 and 1 .2 times the initial guess values (mostly 375 
taken from Paper I), resulting in a significantly lower number of 
misidentified stars. The algorithm automatically identified the 
central radial mode (and therefore provides a mode degree iden- 
tification for all other modes) and determined the central ( + 1 
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radial order around Vmax) large and small frequency separations 
for about 80 % of the stars in our sample. For the remaining 
stars, we iteratively adjusted the parameter ranges by hand to 
achieve a correct identification and finally checked each single 
fit for its plausibility. We determined the most probable param- 
eters and their Icr uncertainties from the marginal distributions 
of the probability density delivered by MultiNest. The large and 
small frequency separations were determined to within 0.4 and 
4.5 %, respectively, for 50 % of our sample and to better than 
1 and 10%, respectively, for 95 % of the stars. The frequency 
of the central radial mode was measured to within 16nHz for 
50 % of the stars and better than twice the frequency resolution 
(~39 nHz) for almost the entire sample. 

3.3. The phase shift of the central radial mode 

For the remaining analysis we consider the asymptotic relation 
for the three central radial orders alone. For radial modes, the 
relation is 



Vcfi 



Av, (n + e'J. 



(5) 



Knowing Vc,o and Avc we can easily determine the phase shift 
of this mode as e/, = (vcfl/Avc) modulo 1. For the Sun, e' is 
small (i.e. closer to than to 1) and per convention e is defined 
as e' -t- 1 to obtain th e correct radial order of the solar p modes. 
IWhite et al] (1201 Ibh showed that e changes smoothly from the 
Sun to the red giants. Hence, to obtain values of e for the red 
giants that are consistent with the solar value we define. 



e - 



e' + 1 if e' < 0.5 and Av > 3 /jHz 
e' otherwise 



(6) 



The top panel of Fig.|4]illustrates the phase shifts of the cen- 
tral radial modes for our sample of 923 red giants as a function 
of Avc (note the logarithmic abscissa). Error bars indicate the 
uncertainties determined from the uncertainties in Avc and Vcfi. 

405 The stars clearly divide into different groups. We subdivided the 
stars to fall above and below a certain limit (dotted line), which 
was initially given by a linear fit (in the e^ - log Av^ plane) to a 
number of stars from the upper population and arbitrarily sh ifted 
by -0.05 in ec. This works fine for stars with Avc above about 

410 3 juHz. For the stars below 3 yuHz, the larger uncertainties blend 
the different populations and we only assigned those to the upper 
group that are within 2cr of the initial (unshifted) fit. The stars 
above the limit presumably belong to the ascending giant branch 
and are indicated by red circles. The stars below the limit are 

415 presumably core He-burning stars, which we further subdivide 
into the following groups. The red clump stars (blue triangles) 
are expected to cover only a narrow range in Avc. They are low- 
mass stars that suffered from electron degeneracy in their RGB 
phase and ignited helium burning in a flash once the core mass 

420 had reached a critical limit. This common core mass means that 
they all produce about the same amount of energy, which ex- 
plains why they span a very narrow range of luminosities and 
radii (and due to a similar mass also Avc). Stars with A vc above 
about 5 juHz (green squares) are secondary-clump stars (iGirardil 

425 1 19991) that are too massive to have undergone a helium flash and 
ignite He in a non-degenerative way, which allows for a range 
of core masses and therefore of luminosities and radii (and Avc). 
There is no counterpart to the secondary-clump stars at the low- 
mass end of the RC below about 3.2 //Hz (yellow diamonds) be- 

430 cause stars with an even lower mass than RC stars are not yet 
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Fig. 4. Various seismic parameters as a function of Av^. Top 
panel: The phase shift e^ of the central radial mode. The stars 
divide into two clear groups, with red circles indicating H- shell 
burning stars. The remainder are core He-burning stars divided 
into stars that presumably belong to the red clump (blue trian- 
gles), the secondary clump (green squares), or the asymptotic 
giant branch (yellow diamonds). The dotted line is an arbitrary 
limit to divide the different groups. Middle panel: The average 
period spacing of mixed dipole modes with the classification 
from the top panel. The dashed line indicates the minimum sep- 
aration that can be measured between two peaks around v^o in 
a 600d-long time series. Bottom panel: Ratio between central 
small and large frequency separation. Colours are only given in 
the online version. 



old enough to having reached this late stage of red-giant evolu- 
tion. We identify them as stars that presumably have akeady left 
the RC to climb the AGB. However, e^ is quite uncertain in this 
domain and some of the AGB stars are likely misidentified. 

To check if our distinction of the evolutionary stages based 
on the phase shift of the central radial mode does actually 
make sense, we follo wed the idea of Bedding et al. (201 1) and 
iMosser et al.l (1201 lah of determining the average period spacing 
of the mixed dipole modes. To do so, we used our mode identifi- 
cation as a template to cut out +2 radial orders (around Vc,o) of / 
- and 2 modes from the power density spectra. The remaining 
signal should then only include dipole modes (and / = 3 modes 
in some cases). We then split the spectrum (now converted to 
periods) into four consecutive parts (two above and two below 
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445 Vcfl), each convering the frequency range between a radial mode 
and the next / - 2 mode, and computed the highest peak in the 
autocorrelation function (above twice the frequency resolution) 
of each part, provided there were at least two well-separated 
peaks present exceeding nine times the background level. The 

450 final average period spacing and its uncertainty were determined 
from the individual measurements and their standard deviation, 
respectively. We were able to measure average period spacings 
for about 2/3 of the stars in our sample. They are given in the 
middle panel of Fig. HI colour-c oded according to the classifi- 

455 cation based on the phase shift. IChristensen-DalsgaardI (|20 1 lb 
argued that He burning causes a convective core where the g 
modes are excluded, resulting in a period spacing in RC stars 
th at is greater than thos e of their RGB counterparts. According 
to'Bedding et al.l (l2011h . the p-dominated dipole modes are sep- 

460 arated by about 50 and 200 seconds in RGB and RC stars, re- 
spectively. The secondary clump stars cover roughly the range 
in between, but at higher Avc than the RC stars. We see that 
our identification of the evolutionary stage based on e^ is cor- 
rect for most of the stars in the sample. About half of the stars 

465 whose classifications do not coincide are located in the region 
where it is difficult to separate RGB from RC stars in the e^ - 
Avc diagram (around the dotted line in Fig.|4). For the remain- 
der, it is likely that we measured multiples of the actual period 
spacing because the individual dipole modes were only poorly 

470 resolved. For stars with Avc below about 3 /iHz the observations 
are not yet long enough to resolve the individual dipole modes 
and significantly longer observations are needed to follow stellar 
evolution towards the tip of the RGB. This clearly shows that the 
phase shift of the central radial mode is as reliable as the aver- 

475 age period spacing of p-dominated dipole modes for determine 
the evolutionary stage of red giants with Avc > 3 //Hz and can 
even provide at least a statistical criterion for stars with Avc < 
3 juHz to distinguish between RGB and AGB stars, which is not 
yet possible with the period spacing. 



480 3.4. The central small frequency separation 



We have also measured the average frequency separation be- 
tween the three radial and / = 2 modes around Vcfi. We show the 
ratio 5vo2/Avc in the bottom panel of Fig.|4] The small spacing 
measures the frequency difference of two modes that have nearly 

485 identical eigenfunction shapes in the outer layers of the star and 
is therefore sensitive to anything that would change the slope of 
the sound speed in the deep interior For stars on the main se- 
quences this includes the slow synthesis of a He core and 6vo2 
can be used as an age indicator It, however, fails as an indicator 

490 of age when the core is not undergoing significant evolutionary 
changes, which occurs during the sub-giant phase where the core 
is still relatively large. On the giant branch, the core is so cen- 
trally concentrated that 6vo2 can once again be used as a measure 
of evolution as it can be sensitive to the state of the H burning 

495 shell (e.g.. White et al. 201 lb). The cores of RGB and RC stars, 
being slightly different, will be affected differently. And this is in 
fact what we see. During RGB evolution, 6vo2 increases contin- 
uously to come back to a value after ignition of core He burning 
in the RC that is on average about 10 % larger than those of their 

500 RGB counterparts. However, different physical phenomena oc- 
curring during different evolutionary stages have similar effects 
on 5vo2 and no clear distinction between them can be made. 




Av, (nHz) 

Fig. 5. Top panel: The phase shift e^ of the central radial mode 
as d erived from th e average Ava with the e(Av) function from 
Mos ser et al.l (l2011 b) superimposed (dashed line shifted by - 
0.05 for better visibility). The classification is taken from Fig.|4] 
with red circles indicating H-shell burning stars and blue trian- 
gles core He-burning stars (including secondary clump and AGB 
stars). Bottom panel: Relative difference between central and av- 
erage Av. Uncertainties are typically well below 1 %. 



3.5. Average vs. central large frequency separation 

To investigate the differences between Fig.[T]and|2] we computed 
the phase shift f or the central radial m ode with respect to the av- 505 
erage Av^ from iMosser et al.l (1201 Ibl) and plot them in the top 
panel of Fig.|5] Obviously, the information on the evolutionary 
stage is hidden in the e^ - Ava diagram indicating that the evo- 
lutionary effect on the radial mode pattern is a local phenomena 
encoded in the fine structure of the m ode sequence. To check if 510 
the approach of lMosser et al.l (1201 Ibh somehow artificially tunes 
Ava so that fa is forced to fall along a perfect relation we re- 
computed our central Av for a broader range of modes and find 
that this is clearly not the case. Now fitting five instead of three 
radial orders, our central Avc is still not a global value, but the 515 
effect of stellar evolution on ec fades with the different popula- 
tions in Fig.|4]blending into each other Interestingly, the differ- 
ence between the average and central Av is only significant for 
core-He burning stars as shown in the bottom panel of Fig.|5] 
Whereas there is no or only a marginal difference for RGB stars, 520 
core He-burning stars have a Ava that is on average about 1.5% 
smaller than their Avc. We note that this opens up the possibil- 
ity to identify the evolutionary stage of a red giant without the 
need to determine a phase shift by comparing the different large 
frequency separations. 525 

3.6. Radial mode curvature of RGB and RC stars 

We find that the difference between the average and central Av 
originates from different shapes of the mode sequences. When 
plotting individual mode frequencies of an RGB and RC star 
with similar Avc in a normalised echelle diagram (top panels of 530 
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Fig-HI, we find that the radial modes follow differently shaped 
sequences. Whereas the radial mode sequence of the RGB star 
(left panel) is "C" shaped (i.e. approximately symmetric around 
the central radial mode), the mode sequence of the RC stars 
(right panel) is more "S" shaped (i.e., antisymmetric around the 
central radial mode). In the case of the symmetric mode se- 
quence, Av determined from the three central radial modes gives 
a very similar value than the average of all modes. Not so for the 
antisymmetric sequences, where we find a significant difference 
between the average and central Av. 

To test if this is a common phenomenon for the different red- 
giant populations, we need to extract individual mode frequen- 
cies for a larger sample of stars. Reliable peakbagging for a large 
number of red giants is a complex task that is not yet solved sat- 
isfactorily and is also beyond the scope of this paper However, 
since we are only interested in an average mode sequence for an 
ensemble of stars, we adopted a very fast and simple approach. 
Using the template 



v„,o 



Vc.O 



+ Av. 



("■^H 



(7) 



with a =^ 0.008 (iMosseretal .1 12011b) and n' being the radial 
order relative to the central radial mode, we could automatically 
estimate the positions of the radial modes. The fit to determine 
the various frequency separations also provides an estimate for 
the average linewidth of Z = and 2 modes, which we find it 
to roughly scale as F ^ 0.026 + 0.12 log Av (with F and Av 
in //Hz). Thi s is in good agreement with what was found by 
ICorsaro et al.l (|2012|) . We then determined the individual mode 
frequencies as the centre of area +3r around the initial guess, 
provided there was a peak present that exceeded nine times the 
background level. 

The result is shown in the bottom panel of Fig.|6] where we 
plot the radial mode frequencies of 48 RGB and 89 RC stars, 
with Avc ranging from 4 to 5 fiHz, in a normalised echelle dia- 
gram. To account for the varying e^ in this range, the abscissa is 
corrected for the phase shift of the central radial mode (as de- 
termined in Sect. 13. 3b . The scatter at the central mode reflects 
the uncertainties of our peakbagging approach and is, with typ- 
ically twice the frequency resolution, significantly larger than if 
we would have fitted proper Lorentzian profiles. Obviously, the 
mode curvature is not only determined by the evolutionary stage 
of the star but the general trend is clear RGB stars have a ra- 
dial mode sequence that is approximately symmetric around the 
central mode and RC stars show a strongly asymmetric or even 
antisymmetric mode sequence. The slope of a linear fit to the av- 
erage mode sequences gives the difference between an average 
(comprising all modes) and central Avc, with a vertical line in- 
dicating no difference. Whereas the difference is negligible for 
RGB stars, the RC stars show a significant difference, with the 
average value being on average about 1.5 % smaller than Avc. 
This rather small difference translates into a significant change 
in fc of typically 0. 14 at Avc = 4.5 //Hz, which is enough to move 
a RC star into the RGB population in the e^ - Ava diagram. 

4. Red-giant models 

To compare the different choices of Av to more physical quanti- 
ties and see whether we can reproduce the observed evolution- 
ary effects on the local fine structure of radial modes, we con- 
structed a sequence of red-giant models. They were computed 
with t he Yale SteUar Evolu tion Code (YREC; Guenther et al. 
119921: iDemarque et al.l 120081) for a near solar composition and 
calibrated mixing length parameter (Z = 0.02, Y - 0.27, aMcr 
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Fig. 6. Top panels: Normalised echelle diagrams of an RGB 
(left) and RC (right) star with 1 = 0, 1, and 2 modes indicated 
by red circles, green triangles, and blue squares, respectively. 
Bottom panel: Radial mode sequence of 48 RGB (open light- 
red circles) and 89 RC (open light-blue squares) stars with Avc 
ranging from 4 to 5 //Hz. The abscissa is corrected for the phase 
shift of the central radial mode. Average values per radial order 
are given by blue-filled squares (RC stars) and red-filled circles 
(RGB stars) connected by broad dashed and dash-dotted lines, 
respectively. Error bars indicate the rms scatter Thin lines give 
the slope of a linear fit to the mode sequences, which differ by 
about 1.5%. 



= 1.8) assuming standard solar mixture (iGrevesse et al.llT996l) . 
More details about the c onstitutive physics may be found in, e.g., 
iKaUinger et al.l ( 1201 Oah and references therein. Like most other 
stellar evolution codes, YREC is not able to evolve models with a 585 
mass below about 2.3 Mq through the helium flash. We therefore 
had to restrict our analysis to low-mass models on the ascending 
giant branch and we only evolved a 2.5 Mq model beyond the 
ignition of He burning until the exhaustion of He in the core 
in the secondary clump. The giant branch part of the evolution- 590 
ary tracks are illustrated in the bottom left panel of Fig.|8] The 
adiabatic pulsation spectrum of every model was calculated us- 
ing Guenther' s n onradial nonadiabatic stellar pulsation program 
(lGuenflieill994 . 

4.1. Average vs. central Av in red-giant modeis 595 

Ideally, the large frequency separation should be measured from 
the highest radial order modes, which follow the asymptotic 
spacing best. These modes, however, are usually not accessi- 
ble or are uncertain due to their small amplitudes. An impor- 
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Fig. 7. Top panels: Ratio of the central (left) and average (right) 
Av to the asymptotic frequency separation A for a sequence of 
red-giant models. Dashed lines indicate unity. Bottom panels: 
Same as top panels but in reference to the square root of the mean 
model density, where we adopt a Av© of 134.88/zHz (Paper I). 



600 tant question, therefore, is how well the various measurements 
of Av compare to the asymptotic spacing and to what extent they 
can be used in the asymptotic relation. Unfortunately, we cannot 
measure the sound travel time across a real star, and therefore 
the asymptotic spacing, directly but we can at least ask how ac- 

605 curately the measured Av follows it for stellar models. 

For a reasonable comparison it is important to extract the 
average and central Av from the models in a si milar way as 
from th e observations. We follow the approach of IWhite et al.l 
( 1201 Ibl) . where only the modes within a specific range around 

61 Vmax were used t o measure Av, with Vmax dete rmined from a scal- 
ing relation (e.g.. lKieldsen & Bedding|1995h . We could compute 
the average frequency separation of the modes within the spe- 
cific range, but this would cause unphysical jumps as modes at 
the top and bottom fall in and out of this range when Vmax varies 

61 5 from one model to the next. A more appropriate approach is to 
compute a weighted average of all modes, w ith the weights give n 
by a Gaussian centred on v^ax (as done bv lWhite et al.ll201 Ibl) . 
The frequency range, and therefore the number of modes for 
which Av is determined, can be controlled by the width of the 

620 Gaussian. For the average Av we chose FWHM - 0.74vinax'^'^^, 
which is a good es timate for the actua l width of the observed 
power excess ("e.g.. iMosser et al.]l2010h . Things are more diffi- 
cult for the central Av, since only a very narrow range has to 
be considered. We found reasonable results with the FWHM be- 

625 ing equal to twice the Av scaled from the model mean density, 
essentially giving weight to only the three central modes. 

The ratios of the central and average Av to the corresponding 
asymptotic spacing (as determined from the models' acoustic ra- 
dius) are illustrated in the top panels of Fig.|7]for models ranging 

630 from 1.0 to 2.5 Mq. We see that both central and average Av un- 
derestimate the asymptotic spacing by typically 1 to 3 %, but 
interestingly there is a clear trend to approach A with increas- 
ing mass. This is because the mode sequence of low-mass stars 
is more curved than those of high-mass stars, i.e. the spacing 

635 of the "unobserved" highest radial order modes differs stronger 
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Fig. 8. Top panel: Phase shift of the central radial mode for a 
sequence of 1.0, 1.5, 2.0, and 2.5 M© ascending RGB models. 
Descending RGB and core He-burning models are shown only 
for the 2.5 Mq track. Bottom left panel: The models from the 
top panel shown in a theoretical HR-diagram with the models 
shown in the bottom right panel indicated with black symbols. 
Bottom right panel: Normalised echelle diagram showing the ra- 
dial modes of 1.0 (red diamonds) and 2.5 Mq (blue circles) RGB 
models and 2.5 M0(grey squares) RC models. The blue and grey 
lines indicate a linear fit to the RGB and RC mode sequences, 
respectively. The grey and red symbols are vertically shifted for 
better visibility. Colours are only given in the online version. 



from a central or average Av (see below). Even though the cen- 
tral Av is more sensitive to mass, i.e. it covers a larger range of 
the ratio Av/ A, it is almost always closer to A than is the average 
Av. This is especially true for the red-giant regime. If we go to 
less evolved sub-giant and main-sequence stars with Av above 
20;uHz (which are not shown in Fig.|7]i, the central and average 
Av tracks converge and eventually fall on top of each other One 
might expect that the deviation between "observed" and asymp- 
totic Av is at least partly due to the so-called near-surface effect, 
which is associ ated with deficiencies in m odelling the upper stel- 
lar layers ("e.g.. lGruberbaueret al.ll2012 ^, and references therein). 
But in fact taking this into account would add another ~0.8 % de- 
viation to the ratio Av/A (at least if we assume the near-surface 
effect for red giants to be solar-like in amount and direction). 

The situation is slightly different when we compare the cen- 
tral and average Av with the square root of the mean model den- 
sity in the bottom panels of Fig.|7] The ratio Av/p'^- (in solar 
units) compares well to Av/A for the low-mass models, but due 
to an even more pronounced mass dependency, it approaches 
unity faster and even overestimates the reference value for high- 
mass models. Again, the central Av appears to be a slightly better 
proxy for the Av scaled from the model density than the aver- 
age Av (except for high-mass stars, which comprise only a very 
small fraction of our red-giant sample). We note that our find- 
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660 ings are fully consistent with those of IWhite et alJ (l2011bl) . who 
performed a similar comparison but using ASTEC models with 
somewhat different input physics. 

Even though we cannot investigate the difference between 
the RGB and RC population for low-mass models, the secondary 

665 clump models do not show any significant difference from their 
RGB counterparts in the comparison with the asymptotic spac- 
ing. Only in the comparison with the model density do we find 
the secondary clump models to differ from the RGB models, 
indicating that the relation between model density and large 

670 frequency separation (i.e., the Av scaling relation) is different 
for RGB and core He-burning stars, as already suggested by 
lMighoetan(l2012h . 

4.2. The phase shift and mode curvature in red-giant modeis 

Finally, we tested whether our models can reproduce the ob- 

675 served evolution of the phase shift and the radial mode curva- 
ture. The top panel of Fig. [8] illustrates the phase shift e^ of the 
models' central radial mode, determined in exactly the same way 
as from the observations. Apart from a systematic shift towards 
lower Ec (which is likely due to the near-surface effect), the RGB 

680 models follow the same linear trend with Avc as the observa- 
tions, with negligible influence from mass. More interesting here 
is the "evolution" of e beyond the ignition of core He burning 
of the 2.5 Mq models. At the tip of the giant branch the stellar 
envelopes start to contract and e^ increases again with increas- 

685 ing Avc but at a slightly lower rate than the RGB models were 
previously decreasing. When the models settle in the secondary 
clump at Avc around 9 fjHz, e^ is about 0.05 smaller than those of 
their RGB counterparts. This is about the same difference we ob- 
serve for many red giants in the corresponding Avc range. Note 

690 that this simple monotonic behaviour of fc is not found for less 
evolved sub-giant and main-seq uence stars, where stellar mass 
has a large influence on e (e.g.. lWhite et al.l[2011bh . This is al- 
ready the case for the 2.5 Mq models with Avc above ~10juHz. 
However, they evolve very fast and are therefore quite rare so 

695 that we do not find a single one in our sample. 

The different mode curvatures are also quantitatively repro- 
duced by the models. This is shown in the bottom right panel 
of Fig. [8] where we plot the radial modes of RGB (blue cir- 
cles) and secondary clump (grey squares) for 2.5 Mq models 

700 with Av ranging from 8 to 9 juHz (indicated by black symbols in 
the HR-diagram). For comparison we also plot the radial modes 
of 1.0 Mq RGB models. Even though these models do not rep- 
resent our observed sample best (for that we would have to in- 
clude 1 .0 Mq RC models, which we cannot compute) the general 

705 trend is clear The mode curvature is largely defined by the model 
mass, although the evolutionary stage has some influence, shift- 
ing the lower order modes towards higher frequencies and the 
higher order modes towards lower frequencies when changing 
from the RGB to the secondary clump. The effect is small but 

71 results in a difference between central and average Av of about 
0.5 % and is therefore accessible in the observations. 

5. The influence of the stellar mass 

We have demonstrated that the various parameters we have cho- 
sen to describe the acoustic spectrum of red giants are affected 
71 5 by the evolution of the stellar core, but it can be expected that 
they also depend on other properties of the star. One of these 
is the stellar mass, which we can estimate via seismic measure- 
ments. Following the approach of Paper I, we compared the mea- 
sured Vmax and Av to those of an extensive grid of stellar models. 
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Fig. 9. Reduced phase shift and central small frequency separa- 
tion for RGB (circles) and core He -burning stars (triangles) with 
the symbol colours coded according to stellar mass. Inserts show 
the mass dependency of the various parameters along with a lin- 
ear fit (red lines). In the insert of the top panel, RGB and RC 
stars are indicated by grey circles and black triangles, respec- 
tively. Note that stars with Avc below 1 .5 fiHz are excluded here 
because their uncertainties are too large. Colours are given only 
in the online version. 



where the seismic parameters of the models were determined 
from scaling relations. A problem in the original approach was 
the ambiguous evolutionary stage of the stars, leading to a rela- 
tively uncertain position in the HR diagram in the region that is 
populated by both RGB and RC stars. Knowing the evolutionary 
stage of the individual stars, we can disentangle the contribu- 
tions of RGB and RC models to the fit. This results in better 
constrained effective temperatures and luminosities but does not 
greatly affect the determination of the stellar mass (or radius). 
With the present observations, the mass is typically constrained 
to within about 5 % assuming that the seismic scaling relations 
are correct. This is not e ntirely true as can be seen from Fig.|7](or 
e.g.. Mig lio et al.ll2012l) and an error of about 2 % in the Av scal- 
ing relation itself adds another about 4 % uncertainty, although 
a detailed error analysis is not of particular interest here and be- 
yond the scope of this paper. 

In the top panel of Fig.|9]we show the influence of the stel- 
lar mass on the phase shift. To correct for the relation between 
Cr and Avc, we fitted a linear function to the RGB stars in the 
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fc - Avc diagram and define fred = ec/(0.58 + 0.66 log Avc). 

740 Our fi t is in good agreement with that found by Corsa ro et al.l 
( 1201 2|) in their study of cluster red giants in the Kepler field 
but slightly different from the relation adopted by Mosser et al.l 
(,2011b.) based on the shorter series from CoRoT. For RGB stars 
we find no significant correlation between Erei a nd M (see the 

745 grey symbols in the inset), which is consistent with' Corsaro et alJ 
(2012 ) and to the study by W hite et al. (2011a) in their extensive 
analysis of stellar model frequencies. The core He-burning stars 
show a slight correlation with M, which does not fully explain 
the dispersion in fred- The dispersion is significantly larger for 

750 RC stars (black symbols) than for RGB stars and can be as large 
as 30 % of Avc. This also explains why we do not find clear / - 
and 2 ridges in the right panel of Fig.|2] The explanation for 
this is unknown but chemical composition might be a good can- 
didate. More interesting here is the gradient of increasing mass 

755 with increasing Avc, which results from the fact that all RC stars 
have a very similar radius (around 1 1 Rq, see Paper I). 

The dependence of 5vo2 on the stellar mass is expected and 
has already been reported for red giants qualitatively from ob- 
servations (Huberetal. 2010) and from theoretical considera- 

760 tions (M ontalban et al.l 1201 0) but has not yet been quantified 
from actual measur ements (apart from the very recent study by 
ICorsaro et al.ll2012l) . From Fig. we see that 6vq2 does not scale 
linearly with Avc, which is presumably due the evolution of the 
burning H shell. To account for this we fitted a power law to the 

765 RGB stars and found that a scaling with Avc **'''* removes the evo- 
lutionary effect best. In the middle and bottom panels of Fig.|9] 
we plot the normalised quantity 6vq2I^vJ^'^'^ for RGB stars and 
core He-burning stars, respectively, with the symbol colour indi- 
cating the stellar mass. Both populations show a clear trend with 

770 5vo2/Avc'^'^'* decreasing for increasing mass. The dependency is 
quantified by a linear fiia + b M (see red lines in the insets), with 
the coefficients a and b being equal to 0.188 and -0.029 for the 
RGB stars and 0.210 and -0.034 for the RC stars, respectively, 
and M in solar units. We note tha t the coefficients pres ented here 

775 differ from those determined by ICorsaro et al.l (l2012l) . which is 
because they assumed a linear relation between the large and 
small separation. However, if we compare our results with theirs 
we derive estimates for (Jv()2 that agree within a few per cent. 

6. Summary and conclusions 

780 When red giants evolve from stars burning H in a shell to stars 
that also bum He in their core, they largely reconfigure their 
interior structure but still populate the same region in the HR- 
diagram and are therefore difficult to distinguish. However, re- 
cent studies have shown that the mixed p/g dipole modes probe 

785 the stellar core and that their period spacings give a particularly 
clean separation of the two groups. In this paper we investigated 
the radial p-mode spectrum of red giants and how it is affected 
by the structural changes due to stellar evolution. Even though 
the radial modes do not probe the core itself, the p-mode cavity 

790 reaches deep into the star and at least its lower boundary should 
be sensitive to the evolution of the stellar core. 

We have analysed high-precision photometric time series 
from the first 600 days of Kepler observations for a large sam- 
ple of red giants. We have measured the local large frequency 

795 separation + 1 radial order around the centre of the power excess 
and determined the phase shift of the central radial mode from a 
local development of the asymptotic relation. Plotting this local 
pair of seismic observables (Avc, ec) against each other reveals 
different populations of stars, with the RC and secondary clump 

800 stars having a systematically smaller fc than their counterparts 



on the RGB. There are also a few stars at low Avc which do not 
fit into the picture of RGB stars, which we identify as stars that 
have already left the RC to climb the AGB. A crosscheck with 
the mixed mode period spacings shows that our classification in- 
ferred only from the three central radial modes is as reliable as 805 
the mixed mode period spacing to distinguish between red giants 
in different evolutionary stages. 

We found a tight correlation between ec and Avc for RGB 
stars. It appears that the difference in e^ between the RGB and 
clump stars becomes smaller and eventually indistinguishable 810 
if we use an average of several radial orders, instead of a lo- 
cal, i.e. only around the central radial modes, large separation 
to determine e^. With this phenomena apparently being caused 
by different large separations, the average Ava turned out to be 
systematically smaller than the central Avc in RC stars whereas 81 5 
there is no significant difference for the RGB stars. Clearly, the 
information on the evolutionary stage is encoded locally, more 
precisely in the shape of the radial mode sequence. This turns 
out to be, on average, approximately symmetric around the cen- 
tral radial mode for RGB stars but asymmetric for RC stars. In 820 
the case of the symmetric sequence, one will always measure 
about the same Av as long as roughly the same number of modes 
are considered above and below the central radial mode. This is 
not so for the asymmetric sequence of RC stars, which results in 
different estimates of Av (that differ on average by about 1.5 %) 825 
depending on whether all modes or only the central modes are 
taken into account. It is also worth mentioning that neither an 
average Av^ nor the central Avc used in the present analysis is 
an ideal parameter since none of them is an exact representation 
of the asymptotic large spacing. However, according to stellar 830 
models, Avc is an at least as good and mostly better proxy for 
the asymptotic spacing than Ava and should therefore be pre- 
ferred in an asymptotic analysis of the observed frequencies. We 
computed radial mode frequencies for a sequence of red-giant 
models and found them to qualitatively confirm our findings for 835 
both the evolution of the phase shift of the central radial mode, 
as well as the with stellar evolution changing shape of the radial 
mode sequences. We also found a clear signature of the evolu- 
tionary stage in 6vq2 and quantified the mass dependency of this 
seismic parameter 840 

Summarising, we found a reliable and relatively easy-to- 
access indicator for the evolutionary stage of red giants based 
only on the three central radial modes. This is of particu- 
lar interest for red giants with very weak dipole modes (e.g., 
JMosser et al.ll2012a) or bright red giants, which are not accessi- 845 
ble to Kepler (or CoRoT) and for which we cannot resolve the 
mixed mode spectrum in the near future and therefore measure 
their period spacing. One of them is e Oph, which is a bright and 
well studie d red giant (e .g., Kallinger et al. 2008). After the ini- 
tial MOST dWalker et al . 2003) o bservations in 2005, e Oph has 850 
been re-observed twice (Ka llinger et al.ll201 ll) so that the data 
are now sufficient to accurately measure the phase shift of the 
central radial mode and to identify the star as a rather massive 
secondary clump star (Kallinger et al. in preparation), which is 
an important input for subsequent modelling. 855 

Acknowledgements. The authors gratefully acknowledge the Kepler Science 
Team and all those who have contributed to making the Kepler mission possi- 
ble. Funding for the Kepler Discovery mission is provided by NASA's Science 
Mission Directorate. TK and JDR are supported by the FWO-Flanders un- 
der project O6260 - G. 0728. 11. SH acknowledges financial support from the 860 
Netherlands Organisation for Scientific Research (NWO) 

References 

Aerts, C, Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, 



11 



Kallinger et al.: Evolutionary influences on the structure of red-giant p-mode spectra 



ed. Aerts, C, Christensen-Dalsgaard, J., & Kurtz, D. W. 

865 Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA Special Publication, 
Vol. 1306, ESA Special Publication, ed. M. Fridlund, A. Baglin, J. Lochard, 
& L. Conroy, 33 
Beck, P. G., Bedding, T. R., Mosser, B., et al. 2011, Science, 332, 205 
Beck, R G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 

870 Bedding, T. R. 201 1, ArXiv e-prints (1 107.1723) 

Bedding, T. R., Huber, D., Stello, D., et al. 2010, ApJL, 713, L176 
Bedding, T. R., Mosser, B., Huber, D., et al. 201 1, Nature, 471, 608 
Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 
Brown, T. M. & Gilliland, R. L. 1994, ARA&A, 32, 37 

875 Chaplin, W. J., Elsworth, Y, Isaak, G. R., Miller, B. A., & New, R. 1999, 
MNRAS, 308, 424 
Christensen-Dalsgaard, J. 2004, Solar Physics, 220, 137 
Christensen-Dalsgaard, J. 2009, ArXiv e-prints (0912.1405) 
Christensen-Dalsgaard, J. 2011, ArXiv e-prints (1106.5946) 

880 Corsaro, E., Stello, D., Huber, D., et al. 2012, ApJ (submitted) 

De Ridder, J., Barban, C, Baudin, P., et al. 2009, Nature, 459, 398 

Demarque, P, Guenther, D. B., Li, L. H., Mazumdar, A., & Straka, C. W. 2008, 

ApSS, 316, 31 
Dupret, M., Belkacem, K., Samadi, R., et al. 2009, A&A, 506, 57 

885 Dziembowski, W. 1977, Acta Astron., 27, 95 

Dziembowski, W. A., Gough, D. O., Houdek, G., & Sienkiewicz, R. 2001, 

MNRAS, 328, 601 
Feroz, P., Hobson, M. P, & Bridges, M. 2009, MNRAS, 398, 1601 
Garcia, R. A., Hekker, S., Stello, D., et al. 201 1, MNRAS, 414, L6 

890 Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010, PASR 122, 
131 
Girardi, L. 1999, MNRAS, 308, 818 

Grec, G., Possat, P., & Pomerantz, M. A. 1983, Sol. Phys., 82, 55 
Grevesse, N., Noels, A., & Sauval, A. J. 1996, in Astronomical Society of the 

895 Pacific Conference Series, Vol. 99, Cosmic Abundances, ed. S. S. Holt & 

G. Sonnebom, 1 17 
Gruberbauer, M., Guenther, D. B., & Kallinger, T. 2012, ApJ (in press); 

arXiv: 1202.2330 
Guendier, D. B. 1994, ApJ, 422, 400 

900 Guendier, D. B., Demarque, P, Buzasi, D., et al. 2000, ApJL, 530, L45 

Guendier, D. B., Demarque, P., Kim, Y.-C, & Pinsonneault, M. H. 1992, ApJ, 

387, 372 
Hekker, S., Basu, S., Elsworth, Y, & Chaplin, W. J. 201 la, MNRAS, 418, LI 19 
Hekker, S., Elsworth, Y, De Ridder, J., et al. 201 lb, A&A, 525, A131 

905 Houdek, G. & Gough, D. O. 2007, MNRAS, 375, 861 

Huber, D., Bedding, T. R., Stello, D., et al. 201 1, ApJ, 743, 143 

Huber, D., Bedding, T. R., Stello, D., et al. 2010, ApJ, 723, 1607 

Iben, Jr., I. 1968, ApJ, 154, 581 

Jenkins, J. M., Caldwell, D. A., Chandrasekaran, H., et al. 2010, ApJL, 713, 

910 L120 

Kallinger, T., Gruberbauer, M., Guenther, D. B., Fossati, L., & Weiss, W. W. 

2010a, A&A, 510, A 106 
Kallinger, T., Guendier, D. B., Matthews, J. M., et al. 2008, A&A, 478, 497 
Kallinger, T., Matthews, J., Guenther, D., et al. 2011, in Proceedings of The 61st 

915 Fujihara seminar, ed. H. Shibahashi, ASP Conf. Series (in press) 

Kallinger, T., Mosser, B., Hekker, S., et al. 2010b, A&A, 522, Al 
Kallinger, T., Weiss, W. W., Barban, C, et al. 2010c, A&A, 509, A77 
Kjeldsen, H. & Bedding, T. R. 1995, A&A, 293, 87 
Koch, D. G., Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79 

920 Mathur, S., Hekker, S., Trampedach, R., et al. 201 1, ApJ, 741, 1 19 
Miglio, A., Brogaard, K., Stello, D., et al. 2012, MNRAS, 419, 2077 
Miglio, A., Montalban, J., Baudin, R, et al. 2009, A&A, 503, L21 
Miglio, A., Montalban, J., Carrier, P., et al. 2010, A&A, 520, L6 
Montalban, J., Miglio, A., Noels, A., Scuflaire, R., & Ventura, P 2010, ApJL, 

925 721,L182 

Mosser, B., Barban, C, Montalban, J., et al. 2011a, A&A, 532, A86 
Mosser, B., Belkacem, K., Goupil, M. J., et al. 201 lb, A&A, 525, L9 
Mosser, B., Belkacem, K., Goupd, M.-J., et al. 2010, A&A, 517, A22 
Mosser, B., Elsworth, Y, Hekker, S., et al. 2012a, A&A, 537, A30 

930 Mosser, B., Goupil, M., Belkacem, K., et al. 2012b, A&A (submitted) 
Roxburgh, I. W. 2010, ApSS, 328, 3 

Roxburgh, I. W. & Vorontsov, S. V. 2000, MNRAS, 317, 141 
Sweigart, A. V. & Gross, P G. 1978, APJSS, 36, 405 
Tassoul, M. 1980, ApJS, 43, 469 

935 Vandakurov, Y V. 1967, AZh, 44, 786 

Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP 115, 1023 
White, T. R., Bedding, T. R., Stello, D., et al. 201 la, ApJ, 742, L3 
White, T. R., Bedding, T. R., Stello, D., et al. 2011b, ApJ, 743, 161 



12 



